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Error Evaluation for Difference Approximations 
to Ordinary Differential Equations 


A recently available report reviews various 
one-step and multistep numerical methods for de- 
riving difference equations to approximate ordinary 
differential equations, and describes a new method 
for evaluating the errors which result from that 
approximation. The method is then applied to a 
nontrivial numerical analysis problem and to an 
investigation of “stiff” equations (equations with 
widely separated eigenvalues). 

The method involves relationships between (1) the 
errors introduced by using finite sampling rates and 
(2) the parameters describing the specific numerical 
method used. The approach treats an analog system 
in terms of isolated modes, and approximates the 
digital algorithm for each mode by its order, p, and 
by certain constants, C, Ao,..*,Ap. The results show 
that the approximation causes a shift of each root, 
given by -Ch^X[’ and also an error associated with 
each forcing function given by 


-ChP^ Ai X'coP-»- 
i = 0 


These results are valid for all small products of 
sampling period, h, and analog system natural 
frequency X, 

The method is primarily intended for use in the 
design and analysis of digital filters and simulators. 
In such applications, errors (primarily due to 


roundoff) are generated during the solution of the 
difference equation, in addition to the approxima- 
tion errors discussed above. These solution errors 
must be evaluated through further analysis. 

Notes: 

1. For related work involving the application of 
finite difference approximation to solving partial 
differential equations, see NASA Tech Brief 
B7M0424. 

2. The following documentation may be obtained 
from: 

National Technical Information Service 
Springfield, Virginia 22151 
Single document price $3.00 
(or microfiche $0.95) 

Reference,: 

NASA CR-103015(N71-16573), 

Difference Equations to Approximate 
Ordinary Differential Equations 

Patent status: 

No patent action is contemplated by NASA. 
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